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Properties of random mixed states of dimension N distributed uniformly with respect to the 
Hilbert-Schmidt measure are investigated. We show that for large N, due to the concentration of 
measure, the trace distance between two random states tends to a fixed number D = 1/4 + I/tt, 
which yields the Helstrom bound on their distinguishability. To arrive at this result we apply 
free random calculus and derive the symmetrized Marchenko-Pastur distribution, which is shown 
to describe numerical data for the model of coupled quantum kicked tops. Asymptotic value for 
the root fidelity between two random states, '/F — 2i/A, can serve as a universal reference value 
for further theoretical and experimental studies. Analogous results for quantum relative entropy 
and Chernoff quantity provide other bounds on the distinguishablity of both states in a multiple 
measurement setup due to the quantum Sanov theorem. We study also mean entropy of coherence 
of random pure and mixed states and entanglement of a generic mixed state of a bi-partite system. 
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I. INTRODUCTION 

Processing of quantum information takes place not in a 
Hilbert space but in a laboratory [T] . However, the stan¬ 
dard formalism of density operators acting on a finite 
dimensional Hilbert space TLn proves to be extremely ef¬ 
ficient in describing physical experiments. It is therefore 
important to investigate properties of the set UIn of den¬ 
sity operators of a given size N, as it forms a scene for 
which screenplays for quantum pieces are written. 

In the one-qubit case the set CI 2 forms the familiar 
Bloch ball, but the geometry of the set Hat C for 

TV > 3 becomes much more complex [2]. Interestingly, 
for large N certain calculations become easier due to 
the measure concentration phenomenon: a slowly vary¬ 
ing function of a random state typically takes values close 
to the mean value [ail]. Extending the approach of OE] 
we focus our attention on the distinguishability between 
generic quantum states. 

Two states p and tr, can be distinguished by a suitable 
experiment with probability one if they are orthogonal, 
i.e. Tipa = 0. The celebrated result of Helstrom |7] 
provides a bound for the probability p of the correct dis¬ 
tinction between arbitrary two quantum states in a single 
experiment, 

P < ^[l + T>Tr(p,cr)], (I) 

where the definition of the trace distance reads 

-DTr(p,o-) = iTr|p-cr|. (2) 
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For two orthogonal states this distance is maximal, 
I^Tr = 1, so p = I as expected. The above bound can be 
achieved by a projective measurement onto positive and 
negative part of the Helstrom matrix, T = p — a. 

In the case of measurement on m copies of both states 
the probability of an error decreases exponentially with 
the number of copies. To characterize the optimal dis¬ 
tinguishability rate one can apply the quantum Stein 
lemma [5] , the Sanov theorem |5] and the relative entropy 
between both states, or the quantum Chernoff bound 

The aim of this work is to analyze properties of generic 
quantum states in high dimensions and their distin¬ 
guishability. For two random states distributed according 
to the flat, Lebesgue measure in the set of quantum 
states of dimension ^ 1, we derive the limiting values 
of fidelity, trace distance, relative entropy and Chernoff 
quantity, which allow us to obtain universal bounds for 
their distingushability. Furtheremore, we study average 
coherence of a random state with respect to a given ba¬ 
sis and average entanglement of random mixed state of 
a bipartite quantum system. 

II. RANDOM QUANTUM STATES AND THEIR 
LEVEL DENSITY 

Consider a Haar random unitary matrix U of dimen¬ 
sion NK acting on a fixed bipartite state \(j)) S 'Hn®'Hk- 
Then the state |'0) = U\(j)) is distributed according to 
the Haar measure on the space of pure states, while its 
partial trace a — Tyk\'4>){4>\ is distributed in the set fljv 
according to the induced measure vk [m. To be specific, 
the integral with respect to the induced measure vk on 
a set of quantum mixed states is equal to the integral of 
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the partial trace over a group of unitary matrices of size 
N X K equipped with the normalized Haar measure, 



f{p)di'K{p) 


f{TTKiUm{00\U^))dpH{U). 


( 3 ) 

In the special case of K = N we get the flat Hilbert 
Schmidt measure, vhs = vn- In particular, setting 
K = N = 2 one generates the flat Lebesgue measure 
on the set H 2 of one qubit mixed states equivalent to 
the Bloch ball. Alternatively, such random states can 
be constructed as p = GCyTiGG^ where G stands 
for a random rectangular matrix of dimension N x K 
with independent complex normal entries. In the case 
where K = N large, the flat measure in the set 
leads asymptotically to the Marchenko-Pastur distribu¬ 
tion MV{x) [12]. In the general case this distribution 
depends on the rectangularity parameter c = K/N. For 
c > 1 the density is continuous. 


MVc{x) = ^ - (1 - Vc)2] [(1 -h v^)2 - a;], (4) 


while for c < 1 the measure also contains a singular com¬ 
ponent and reads (1 — c)5o + MVc{x). 

Here x = KX is a rescaled eigenvalue of p. In the case 
where c = 1 the Ai V measu re p.MV has a continuous 
density, AiV{x) = [yjAjx — l)/27r, supported on [0,4]. 

Consider a random quantum state p generated accord¬ 
ing to the flat measure r'ns in Hjv with eigenvalues 
described by the measure p!'^'^ = ^J2i^NXi{p)- For a 
large N with probability one the measure pN converges 
weakly to Pmv- Thus, for any bounded and continuous 
matrix function g of a random state p, the trace of g{p) 
converges almost surely (a.s.) to the mean value. 


We wish to evaluate the trace distance between two 
random quantum states p, tr sampled from the set 
according to an induced measure v^- To this end we 
need to derive the spectral density of the Helstrom ma¬ 
trix r = p — a, for which we use tools of free probability. 
In particular, we are going to apply the free additive con¬ 
volution of two distributions m, written P 1 WP 2 , which 
describes the asymptotic spectrum of the sum of two ran¬ 
dom matrices Xi and X 2 , each described by the densities 
Pi and P 2 , under the assumption that both variables are 
independent and invariant with respect to unitary trans¬ 
formations, Xi —>■ UXiW. As both random states p and 
a are independent and free, the limiting spectral density 
of their difference is given by the free additive convolution 
of the Marchenko-Pastur distribution and its dilation V 
by factor minus one, SMV = In order 

to derive SAiV, it is convenient to use the R transform 
of the XIV distribution, given by: 

R{G{z)) + —^ = z, (8) 

G[z) 


where G{z) is the Cauchy transform of the distribution 
m- Hence, we get that the i?-transform of the MP dis¬ 
tribution reads Rmv{z) = c/(I — z) so the A-transform 
of the symmetrized distribution is given by the sum. 


Rsmv{z) = Rmv{z) — Rmv{—z) = 


2cz 

I-z' 
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Inverting this transform we obtain the desired 
metrized Marchenko-Pastur (SMP) distribution: 


SMVciy) = 


-1 - 4(c - l)c - -k [Ye(y)]^^^ 
2^/^^Ty[Yc{y)]^/^ 


( 9 ) 

Sym- 

( 10 ) 


Tr( 5 (p)) = nJ g dp^^\t) 

where g{t) = limjv_,.oo Ng{X). 


g{t)dpMv{t), 

( 5 ) 


where 


Y,iy) =(2c - If + 9(c + l)j/2 + sVdy 

X a/(2c — 1)3 — y^ -\-[2 — {c— 10)c]?/2. 


( 11 ) 


III. TRACE DISTANCE AND A SINGLE SHOT 
EXPERIMENT 

To show a simple application of Eq. ([^ consider the 
trace distance of a random state p to the maximally 
mixed state p» = 1/N. For a large dimension N the 
value 


Tt:\p- pf= j \t-l\dpf\t), ( 6 ) 

converges to the integral over the AiV measure, so that 
the trace distance behaves almost surely as 

T>Tr(p,P*) > \ j\t- Ci 0.4135. 

Af->oo 2 J 47r 

( 7 ) 


This measure is supported between y± = 
-f lOc -I- (c -I- -\- 2. If c < 1 the 

measure has a singularity at zero. 

In the special case of c = 1, corresponding to the 
Hilbert-Schmidt measure z^hs in Hjv, we obtain 


SMP[y) 


-1 - 32/2 + (1 + 3p (v/3 + 33p2 - 3p4 + 6j/) ) 

1 

2y/dTTy (^1 -k dy d-dy'j'j 

( 12 ) 


This measure was earlier identified as a free commutator 
between two semicircular distributions m and appeared 
in the literature under the name of the tetilla law [ID¬ 
AS a simple application of the SAiV distribution we 
will sow that the Hilbert-Schmidt distance tends to zero 
as N tends to infinity. 
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FIG. 1: Symmetrized Marchenko-Pastur distribution 
(10) for c = 0.2 (a),0.5 (6), 1.0 (c),4.0 (d) denoted by 
the solid (red) line. Histograms represent numerical 
results for the model of coupled kicked tops (67 ^ 


with N = 100 cumulated out of 100 realizations of the 
system. 


To obtain the required result consider two indepen¬ 
dent random states p and cr, distributed according to the 
Hilbert-Schmidt measure and rewrite Eq. ([^ as: 

Trigip - a)) = N J g (13) 

From this follows: 

lip - o-|l2 = Tr(p -af = (14) 


Proposition 1 . Consider two independent random 
states p and a drawn aecording to the HS measure in 
the space of states of dimension N. For large N the 
probability p 2 of correct distinction between both ran¬ 
dom states in a single measurement is bounded by p 2 < 
5 5 (i + w) — 0.7842, while the probability pi of dis¬ 

tinguishing a and the maximally mixed state = 1 /N 

is bounded by pi < ^ -\- — 0.7067. 


The above results improve our understanding [5] of the 
structure of the set Hjv of mixed states of a large dimen¬ 
sion TV. The HS measure is concentrated in an e neigh¬ 
borhood of the unitary orbit, UpW, where U is unitary 
and p is a random mixed quantum state with spectrum 
distributed according to pmv- The width of the orbit 
decreases as - this will be shown in the next section ~ 
while its diameter is given by the distance between two 
diagonal matrices with opposite order of the eigenval¬ 
ues [I?], d = T>Tr(p^,P'^) = /o'a;sign(a: - M)dpMv{x) — 
0.7875, where M denotes the median, dpMV = 1/2. 
A generic state p is located at the distance a = from 
the center p* = T, its distance to the closest boundary 
state p tends to zero, while the distances to the closest 
and to the most distant pure state, \4>i) and | 02 )j tend 
to unity - see Fig. 

To study the case of a large environment, AT 3> TV, 
we take the limit c —>■ oo of Eq. (10). This distribution, 
rescaled by a factor 1/v^ tends to the Wigner semi¬ 
circle law C{y) supported on [—2,2], as its i?-transform 
converges to Rw{z) = z. The average \y\ over the semi¬ 
circle is finite in analogy to but due to the scaling 

factor the typical trace distance with respect to the 
limiting measure Vca tends to 0 - see Appendix 


Now taking the limit TV —>■ oo, the integral with proba¬ 
bility one has a finite value f t^dpsMv(i) = 2, hence due 
to the prefactor TV“^ we get asymptotically the desired 
result. 


idHs(p,a) = ||p-a||2CsyA_^0. (15) 

An analogous result implying that the average HS dis¬ 
tance Uhs (P) P *) between a random state p and the max¬ 
imally mixed state p* for large dimension TV tends to zero 
was earlier obtained in [4]. Although this distance is not 
monotone [5] it admits a direct operational interpretation 

[IB]- 

The asymptotic value of Trjp — crj = Trjrj will be al¬ 
most surely given by an integral over the SA4V distribu¬ 
tion, which yields the generic trace distance between two 
random states. 


> D = i / \y\dpsMv{y) = 7 + - - 0.5683. 


T'Tr(p,Cr) 

AT—>-00 Z / ^ ‘/i 

(16) 

Formulae ([^ and (16) combined with the Helstrom 
theorem imply one of the key results of this work. 


P 



FIG. 2: Sketch of the set Hat of mixed states for large 
dimension: the measure is concentrated along the 
unitary orbit of a generic state p = p^. Its trace 
distance to the center p* of the set is a ~ 0.41 and to 
another typical state cr is l) ~ 0.57. 


Note that the trace distance asymptotically tends to 
a fixed value (16), while the average HS distance tends 
to zero. This is a consequence of the known fact that if 
the taxi distance between two points in is fixed. 
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their Euclidean distance is typically much smaller, since 
shortcuts across city blocks are allowed. More formally, 
one can show that a random point from the unit sphere 
S'” with respect to the taxi metric Li is located at the 
Euclidean distance L 2 1/v^ form the center of the 
sphere. A similar relation between the trace and the 
Hilbert-Schmidt distance holds in the space of quantum 
states and explains result (15). 

To visualize this, we present here, for comparison, the 
average distances between two random points in a unit 
ball in M" distributed according to the Euclidean mea¬ 
sure. The distances with respect to the norms Li, 
and Loo are provided in Tab. |lj Note that the definition 
of the trace distance contains the prefactor 1/2, so the Li 
distance has to be compared with 2L)xr- For large dimen¬ 
sions the full mass of the n-ball is concentrated close to 
its surface, so the Li distance between two typical points 
diverges, while in the case of the set of mixed states, the 
layer of generic states remains distant from the subset of 
pure states and the average asymptotic distance becomes 
finite and deterministic. 


two states \tp),\(j)) € 'Hn ® 'Hn and the function 


= DT,{p,a), (19) 


where p = Ttb\'4^) {'4’\ and a = "Itb\(1>){4>\ represent re¬ 
duced states. For a fixed state \(j)) the above function 
has the Lipschitz property with the Lipschitz constant 
equal to unity. To show this we write 


<l\\p-p'\\u 






f^iii 


( 20 ) 


where p' = Trsl'C))?!- Next we rely on the inequality of 
Fuchs and van de Graaf mm. 

7^Tr(p,p') < (21) 


where the fidelity F between two states reads 


n 

Li 

La 


1 

2 

3 

2 

3 

2 

3 

2 

~ 1.1528 

^ ~ 0.9054 

4r>7r 

0.8151 

3 

~ 1.15428 

§ ~ 1.0286 

0.8549 





00 

CXi 

V2 

0 


TABLE I: Average distances of points in a unit ball of 
dimension n with respect to the norms Li, L 2 and L^a- 


IV. APPLICATION OF LEVY’S LEMMA 


L(p,a) = (Tr|pi/Vi/"|)2. (22) 

Using an elementary inequality \/l — a:^ < \/2 — 2x for 
X G [0,1], we write 

nr,(p, p') < < II IV’) - 10 II. (23) 

The last inequality follows from the Uhlmann theorem 
m relating fidelity to the maximal overlap of purifica¬ 
tions. This consideration gives us, that for a fixed |^) 
the function fdip), |</)) is 1-Lipschitz. 

Using the Levy’s lemma, for a fixed state \(j)) and \ip) 
distributed uniformly on a set of pure states in "Hjv G 
we obtain 


In this section we present the Levy’s lemma, used to 
strengthen our results and to characterize the conver¬ 
gence rate of the trace distance between two random 
states to the mean value D. Note that the key results 
from the previous section, showing that various distances 
between two generic states converge to their determinis¬ 
tic values, are independent of the Levy’s Lemma. 

Lemma 1 (Levy’s lemma) . Let f : ^ K &e a 

Lipschitz-continuous function with Lipschitz constant rj, 

i.e. 


fix) - fiv) 


< nWx - : 


(17) 


( /(l^)> !<(>))-E/(|i/’),l</)) >e) 

( DtAP. - EDTr(p, cr) > 


P 

= P 
< 2 exp I — 


(24) 


9tt^ 


Let us now consider a sequence of fixed states {|</Ar)}Ar, 
such that the empirical distributions of eigenvalues of the 
reduced states crv tend to the Marchenko-Pastur law. 
For a random state \ip) G "Hat with p = "ItB\' f’){'f’\ 

we define an expectation value with respect to the distri¬ 
bution of I'i/i), 


For a random element of the sphere x G S'™ ^ distributed 
uniformly, for all £ > 0 we get 

Pi\f{x)-Ef{x)\>£)<2exp(^-^^^"^. (18) 

Following [3] we are going to use the Levy’s lemma for a 
set of bipartite pure states of a high dimension. Consider 


Aat = E/(|^/), |(/Af)) = EDTr(/0, ctat). (25) 


Unitary invariance of the trace distance and uni¬ 
tary invariance of the distribution of p implies 
that Ajv -G D = j + ^. The triangle in¬ 

equality gives us the following set of inclusion 


relations a^) — A^r 


-k 


An — D 


< e 


c 
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{p: 

that 


DTr(p, o-at) — D < . For dimension N so large 

An — D\ < e this implies the following bounds, 

(^Dti {p,(Tn)) - D < 

DTrip, ctn)) — An 


P 
> P 
= P 


An — D\ < £ 
DTt{p,aN)) — An <e—|A jv — 

^ N\e-\AN-b\f\ 


(26) 


>1 — 2 exp — 




Since |A^r — Zi| = o(l) we arrive at the following result. 

Proposition 2 . Consider two independent random 
states p and a distributed according to the Hilbert- 
Schmidt measure in the set Hn of quantum states of di¬ 
mension N. For sufficiently large N, when |EI?Tr(p, f) — 
D\ < |, their trace distance is close to the number 
H = I + as the deviations become exponentially rare, 


( £)Tr(p,cr)) 


-D 


> £) < 2 exp — 


N^{£/2f 


Qtt^ 


(27) 


The formal statement above can be interpreted that 
the orbit of generic quantum states containing the entire 
mass of the set i^N, represented in Fig. by an e-strip, 
for large N becomes infinitesimally narrow as e ^ 1/N. 
Since the exponent in Eq. (27) behaves as iV^e^, the 


width £ of the strip decays with the dimension faster 
than for any (5 > 0. 


V. SEVERAL MEASUREMENTS, RELATIVE 
ENTROPY AND CHERNOFE BOUND 

Consider the Chernoff function [TOl HI] of two ran¬ 
dom states, Qs{p,cr) = Trp'*cr^“® and the Chernoff in¬ 
formation Q := min^gjQ 2 ] Qs- Based on the discussion 
in Appendix the Chernoff function can be asymptot¬ 
ically represented as a product of two integrals over the 
Marchenko-Pastur measure. 


Qs > [ PdpMvit) X [b ®d/r(vt'p(0 
AT-i-OO J J 


N—foo 

4r(|-s)r(s + i) 


(28) 


7rr(3 - s)r(s + 2) 

The above expression, involving the Gamma function 
F(a;), takes its minimal value at s = so that the generic 
value of the quantum Chernoff information Q reads 

Q = minf lim Qs) = f ~ 0.7205. (29) 

s VTV-foo / J 

The Kullback-Leibler relative entropy reads S'(pllcr) = 
Trp (logp — log (t) . Applying again the reasoning pre¬ 
sented in Appendix]^ we get: 

J{tlogt - t\ogs)dpMv{t)dpMv{s) = 

(30) 


S{p\\a) 


N—>-oo 


In general the relative entropy is not symmetric, but for 
two generic states, belonging to the same unitary orbit 
one has S{p\\a) = S{a\\p). In the case when one of the 
sates is maximally mixed the symmetry is broken and we 


get ^(pIIp*) 


N —>-oo 


> while 5'(/9 *||p) > 1. Results 


N—*-oo 


(29) and (30), compatible with [5], combined with the 
quantum Sanov [9] and Chernoff bounds [iniiis] imply: 

Proposition 3 . Perform m measurements on a generic 
quantum state p of a large dimension. Probability of ob¬ 
taining results compatible with measurements performed 
on another generic state a scales as exp(—3m/2). In a 
symmetric setup, probability of erroneous discrimination 
between both states behaves as exp(—64m./97r^). 


VI. TRANSMISSION DISTANCE AND BURES 
DISTANCE 

The Jensen-Shannon divergence (JSD) characterizes 
concavity of the entropy. Its quantum analogue (QJSD) 
deals with the von Neumann entropy, S{p) = —Trplogp, 
and is equal to the Holevo quantity y, 

QJSD{p, a) := S [{p + a)/2] - [5(p) + S{a)]/2. (31) 

Level density of the sum of two generic random states, p-l- 
CT, is asymptotically described by the MV 2 distribution. 
As the average entropy for AiVc distribution for c > I 
reads [20], 

{S)c = -J X log xdpMV„ix) =clog c, (32) 
using Eq. ([^ it is easy to show that 
QJSD(p,cr) =yt[logt-log2V]d^^‘^(t) (33) 


+ ^ / t[logt- log A][dAi^(t) -k dp%{t)] - • > 

Z J N^oo 

j t'^ogtdpMvAt) +log2-k J tlogtdpMv{t) = y 

JSD is a symmetric function of both arguments but it 
does not satisfy the triangle inequality. However, its 
square root is known m to yield a metric for classi¬ 
cal states, which leads to the transmission distance [22] . 
Dt{p,cf) = -\/QJSD(p, cr) 1/2. 

Squared overlap of two pure states, \{if\(j))\^, can be 
interpreted in terms of probability. A suitable general¬ 
ization of this notion for mixed states, called fidelity, is 
defined in (221. The average fidelity and root fidelity be¬ 


tween random states for small N was analyzed in [23] . 
To cope with the same problem for large dimensions we 
use the Fuss-Catalan distribution [Ml. 


FC{x) = 


f/2y/i [f/2 (27 -k 3V8I - 12a:) ® - 6^] 


I27r 


xl (27 -k 3V8I - 12a:) ■ 


(34) 
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VII. AVERAGE ENTANGLEMENT OF 
QUANTUM STATES 
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FIG. 3: Dependence of average distance between two 
generic states on the dimension N. Dashed (red) line 
shows the Bures distance and solid (black) line shows 
the trace distance. The horizontal lines mark the 
asymptotic values. 


which is supported in [0, 27/4] and describes the asymp¬ 
totic level density of a product pa of two random density 
matrices. Hence, the average root fidelity can be asymp¬ 
totically expressed as an integral over the Fuss-Catalan 
measure, 

^/F{p,a) = ^/X^{pa) J y/x FC{x)dx = ^. (35) 


This important result, compatible with [23], determines 
the limiting behavior of the Bures distance [5], 


Db{p, a) = ^2{l-yjF{p,a]) ^ \/2(l - 3/4) = ^.(36) 


Observe that the known bounds [19] between hdelity, 
Chernoff information and trace distance, F < Q < y/F < 
\/l- D^, are generically satisfied with a healthy margin, 
0.562 < 0.720 < 0.75 < 0.823. 

Making use of the Marchenko-Pastur distribution we 
calculate the average root fidelity between a random state 
p and the maximally mixed state p*. 


y/F{p,p^) J y/idpMv{t) = ^ ~ 0.8488 (37) 


In a similar way one can treat the related Hellinger 

distance [2], Dnip^a) = \J‘2 — 2 TrpW 2 . Making use of 
the average (29) we find the asymptotic behavior. 


Dh{p, cr) 


/ 2-2 


8 

37r 


~ 0.7476 


(38) 


Dependences of averaged values of selected distances 
between two generic states on the dimensionality shown 
in Fig. [^ are consistent with rumors that form the point 
of view of random matrix theory a dozen is already close 
to infinity. 


Generic pure quantum states are strongly entangled 
[1 E]. Integration over the A4P measure may be 
applied to evaluate average entanglement of a random 
pure state \tp) € 'Hn ®'F.n- Its G-concurrence is de¬ 
fined [2S] as G = iV(detp)^/^, where p = Tj:jy[\'tjj){ip\. 
For large N the concurrence converges almost surely to 
(G) —>■ exp[J \og(x)dpMv{x)] = I/e, in accordance with 

m- 

Results developed in the former sections allow us to 
study the overlap between a given pure bipartite state 
\ip) G fljv <8 Hat and the nearest maximally entangled 
state. This overlap, also called maximal entanglement 
fidelity, is related to another measure A/” of entanglement, 
called negativity. The relation is as follows 




max |('0|^+)P —I 
2 |p+)Gn“f 


(39) 


where denotes the set of maximally entangled states 
and Xi are the Schmidt coefficients of the state \ijj), which 
satisfy = 1- For a partial trace, p = Tri|^/) (■!/|, we 

have 


^ = Try/p, (40) 


Now we consider a random pure state distributed accord¬ 
ing to the unitarily invariant Haar measure. Integral (37) 
implies that 


1 

y/N 




8 

Stt ’ 


(41) 


almost surely when V —)■ oo. The above result gives us 
the asymptotic behaviour of the negativity for pure states 





(42) 


Formally we have an almost sure convergence in the 
asymptotic limit, N ^ oo, 

^A7(|V^)(V^|) ^ i (^1^) « 0.3602 (43) 

Average entanglement of random mixed states can also 
be analyzed m- As shown by Aubrun [55] the level den¬ 
sity of partially transposed random states, written p '^^, 
of a bipartite system obeys asymptotically the shifted 
semi-circular law of Wigner, 

- ;^a/4c- (x- c)2. (44) 
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Therefore, assuming 0 < c < 4, the fraction /v of nega¬ 
tive eigenvalues depends on the measure, 


/at = / 7^ y^4:C-{x- cydx 

Jc-2^ Sttc 

darccos — a/4c — 

47r ’ 


(45) 


and tends to 0 for c —>■ 4 in agreement with [5S]. 
thermore, the average negativity, Af{p) = Trjp'^'^l — 
a random mixed state behaves as 

N = [ \/4c — {x — c)'^dx = 

Jc-2^c 27rc 

SV4c — -I- V— 12c arccos 


Fur- 
1, of 


(46) 


127r 


For the flat measure, c = 1, we get = ^ — 0.1955 

and M = i ~ 0.080 respectively, which implies 

that a generic bipartite state is weakly entangled. The 
above numbers can be compared with the values for the 
maximally entangled state, p+ = which 

read fj^{p+) —>■ ^ and Af{p+) = — 1. 


a random pure state \tjj) in the computational basis, (or 
the basis state |1)(1| in a random basis) we arrive at a 
random projector matrix uj = \ip){'ijj\ = uP'. The coarse 
graining channel describes decoherence, as it sends any 
state p into a diagonal matrix, 


$cg(p) = XI l*)(*l^l*)(*l = diag(p). (48) 


Applying this channel to a projector w we obtain a matrix 
with classical probability vector, qi = wn = |(i|^/>)P, at 
the diagonal. If l"^) represents a complex random vector, 
then the joint probability distribution P(( 7 ) is given by 
the Dirichlet distribution (471 with s = 1, while random 
real vector leads to the case s = 1/2. 

Integrating out TV — 1 variables from the Dirichlet dis¬ 
tribution (47) one obtains the distribution of a single 
component x = qi oi the probability vector. For a fixed 
dimension N these distributions read, for s = 1/2 and 
s = 1, respectively. 


Pl/2{x) 

Pi{x) 


r(f) (1 - a:)(^-3)/2 

r(X^) 

{Nx)^-^. 


(49) 


VIII. CLASSICAL PROBABILITY VECTORS 
AND AVERAGE COHERENCE 


Let q = {qi,...qpf) be a normalized probability vec¬ 
tor, so that qi> 0 and ^ ■ qi = 1. A vector q belongs to 
the probability simplex of dimension N —1, in which one 
defines a family of symmetric Dirichlet measures, param¬ 
eterized by a real index s. 


N 


Ps{q) = Cs X - 1 (9112 ■ • • QnY 


(47) 


Here Cs denotes a suitable normalization constant, which 
is dimension dependent. Note that the flat measure in 
the simplex is recovered for s = 1, while the case s = 1/2 
corresponds to the statistical measure [5]. 

Both of the above measures can be related to ensembles 
of random matrices. The flat measure, s = 1, describes 
distribution of squared absolute values of a column (or 
a row) of a random matrix distributed according to the 
Haar measure on the unitary group, while the statistical 
measure, s = 1/2, corresponds to random orthogonal ma¬ 
trices m- It is known [29] that eigenstates of matrices 
pertaining to circular unitary and circular orthogonal en¬ 
sembles are distributed according to the Haar measure on 
unitary and orthogonal groups, respectively. Therefore, 
these distributions characterize statistics of eigenvectors 
of random unitary matrices pertaining to circular uni¬ 
tary ensemble (CUE) and circular orthogonal ensemble 
(COE), respectively [30] . 

It is also possible to relate these measures with ensem¬ 
bles of random quantum states. Writing a projector on 


In the asymptotic case, iV —> oo, the above distributions 
converge, after an appropriate scaling, to the xt distri¬ 
bution with the number ly of degrees of freedom equal to 
1 (orthogonal case) and 2 (unitary case) |5S]: 


A/2(.)^X?(t)^^e-V2, 

p^t) = 2xim = e-Y 


(50) 


Here, t = Nx, denotes a rescaled variable, such that 
(t) = 1, due to the fact that an element (x) is pro¬ 
portional to 1/7V. The former formula, also known as 
the Porter-Thomas distribution, describes statistics of a 
squared real gaussian variable, while the latter one de¬ 
scribes distribution of the squared modulus of a complex 
gaussian variable. 


Using the distributions (50) we calculate the asymp¬ 
totic mean Li distance between two probability vectors 
p and q, distributed according to the statistical D®, and 
the flat measure D^, respectively. Furthermore, we find 
their root fidelity: 


y/np;q) = Yl VpPb- (51) 

i 


This quantity is also known as the Bhattacharya coeffi¬ 
cient [2]. Finally, we compute the root fidelity and the 
Bures distance from the maximally mixed vector p*. We 
get the following results for the root fidelity \fF^ and 
ppP averaged over statistical and flat measures, respec- 















tively, 


(VF^{p,P*)) ^ JVtPi/2{t)dt = « 0.7979, 

^ IViPiit)dt = ^ « 0.8862, 

\/F®(p, JV^Pi/2it)Pi/2is)dtds = “ ~ 0.6366, 
^F^(p, —>■ JV^Pi{t)Pi{s)dtds = ^ ~ 0.7853. 


(52) 

Observe that for both measures the average root fidelity 
between two random probability vectors asymptotically 
equals the average hdelity with respect to the uniform 
vector, {yjF{p,q)) = {^/F{p,p^)f = (F(p,p*)). 

The above numbers can be directly compared with 
results obtained in Sec. |VI| for the quantum case. For 
instance, the average root fidelity between two random 
classical states distributed according to the flat measure, 
(0.78), is larger than the analogous result, (0.75), for 
quantum states. Hence the average Bures distance satis¬ 
fies the reversed inequality, {DB{p,q)) < {DB{p,cr)). 

Mean Li distances between classical states averaged 
over the probability simplex can be expressed as inte¬ 
grals over the measures (49). For large N one can use 
asymptotic measures (50), which yield the following re¬ 
sults, 

^ {D%p,p*)) J\t- MPi/2it)dt = 

i (£>f(p,p*)) J \'t- l|-Pi Wdt = K 

\{D%P,q)) J |t-s|Pi/2(t)Pi/2(s)dMs = ^, 

\{D^iP,Q)) J\t-s\Piit)Pi{s)dtds = 


(53) 


To make a direct comparison with the results obtained 
in the quantum case easier, each formula above contains 
the same prefactor 1 / 2 , present in the definition ([^ of 
the trace distance. 

In order to quantify coherence of a state p with respect 
to a given basis one uses the relative entropy of coherence 
EH, defined as the difference between the entropy of the 
diagonal of the density matrix and its von Neumann en¬ 
tropy. 


Crel.entip) = S{d>CG{p)) - S{p). 


(54) 


It belongs to [0,logA^] and is equal to the increase of 
the entropy of a state p under the action of the coarse 
graining channel (48). 


The maximal value, Crei.ent = log N, is achieved for 
contradiagonal states [31], for which pa = l/N. If t/min 
denotes the unitary matrix of eigenvectors of p, so that 


Pd = ULinPF^niin is diagonal, then pc = Uf^^^pUraa.^ is 
contradiagonal in the related basis t/max = U-minF for 
any unitary complex Hadamard matrix F, which satisfies 
FF^ = 1 and |Fy | = l/N 


For any pure state the second term in (54) vanishes. 


so for a complex random pure state of a large dimension 
N its average entropy of coherence equals 


(C'rel.ent(IV’c))) = {S{p))cUE = 

= - J tlogtPi{t)dt = log IV - (1 - 7 ), 


(55) 


where 7 « 0.5772 denotes the Euler constant. An analo¬ 
gous result for a random real pure state reads 


(C'rel.ent(IV'M))) = (-S'(p)) COB = 


= - y t log tPi /2 (<)dt = log - (2 - 7 - log 2), 


(56) 


where the average entropies were studied first by 
Jones |34j . Thus a random state is characterized with 
a high degree of coherence, close to the maximal one, 
logA^. Interestingly, this statement holds almost surely 
for any choice of the basis, with respect to which the de¬ 
coherence takes place. We also note that these results 
are compatible with the ones presented in |35j when we 
let N —>■ 00 . 

For comparison let us consider a random mixed state p 
distributed according the Hilbert-Schmidt measure. Its 
average entropy is close to the maximal, as it reads | 2 ] 
for large dimensions, {S{p))hs = ^ogN — 1 / 2 . There¬ 
fore, its relative entropy of coherence, equal to the en¬ 
tropy gain induced by the decoherence channel (48) can¬ 
not be large. In the generic bases the diagonal ele¬ 
ments of a complex Wishart matrix are asymptotically 
distributed according to the Xv distribution with ly = 2N 
degrees of freedom, the entropy of the diagonal behaves 
as S'(diag(p)) ~ logA^— 1/2A^. Therefore, the relative 
entropy of coherence of a random mixed state p asymp¬ 
totically tends to a constant, 

(Crei.e„t(p)) = (5(diag(p))) - {S{p))hs ^ 1 / 2 , (57) 

characteristic to a random mixed states in a contradiago¬ 
nal form, for which Crei.ent is maximal. Thus the average 
coherence of random states decreases with their purity, 
as expected: It behaves asymptotically as logA^ for a 
random pure state and it is equal to a constant for a ran¬ 
dom mixed state. This result is compatible with the one 
presented in Tab. II in [36] . when we multiply the values 
in the second column by log N. 

Note that result (55) can be used to explain recent 
findings of [37]. The mean entropy {Sdia^ij)) of the di¬ 
agonal of a quantum state averaged over a sufficiently 
long time r tends to the average over the ensemble of 
random states, which in the case of complex pure states 
tends to the mean entropy (55) over the flat measure 


on the simplex. It is equal to logA^ — (I — 7 ) and be¬ 
comes larger for mixed states. As the entropy of the 
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time averaged state Sp is limited by logA^, their differ¬ 
ence AS' = Sp — (Sdiag(T)) for a generic unitary evolution 
satisfies the bound AS < 1 — 7 , in agreement with [55] , 
However, taking an initial real pure state pn and restrict¬ 
ing attention to purely imaginary generators, H = iA 
with A real antisymmetric, A = —A^, one obtains an 
orthogonal evolution matrix, O = = e“^, such that 

the state = OprO^ remains real. In such a case the 
average entropy of the diagonal is equal to the entropy 
averaged over the statistical measure (56 1 , so that the 


bound for the entropy difference becomes weaker. 


AS < 2 - 7 -log 2 « 0.7297 


(58) 


Next we wish to calculate the Li coherence of a generic 
quantum state 


Cl, =Y.\P^,A■ ( 59 ) 


In order to achieve this observe first that for any N > \ 
and a random pure state p = \4‘){4>\ distributed according 
to the Haar measure, we have: 

I Pit I = ( 60 ) 

Next, we note that the vector q = p,..., 
pertains the A-dimensional symmetric Dirichlet distri¬ 
bution. For real states we get q ~ Pi/ 2 [q) and for 
complex states we get q ~ Pi{q)- Next, we note that 
®(|Pii)l = calculate the above expecta- 

tional value we will use a formula for mixed moments of 
the Dirichlet distribution. For a random vector q dis¬ 
tributed according to the Dirichlet distribution with pa¬ 
rameters given by the vector a we have |39j 



r(Z)fcl •ai + Pi) 


X 


n 




r(ai -I- Pi) 

T{a,) 


(61) 


To calculate the expected coherence of a pure state we 
set ai = 1 for complex states and at = 1/2 for real 
vectors. In both cases we set /3i = /32 = 1/2 and Pi = 0 
for i > 3. We also note that for large N, we have Cl^ ~ 
N{N — l)(|py I). We obtain 

(« = (A- 1 )^, 

^ (62) 
Cliiic) = {N - 1)^. 


Values obtained above are related to the limiting values of 
the root fidelity for classical probability vectors presented 
in Eqn. (52). This happens because in the asymptotic 
scenario the distribution of the probability components 
and non-diagonal elements of a pure state, differ only by 
a scaling factor. 

These values can be compared with the maximal value 


attained for pure state in a contradiagonal form, p^ for 
which \{pc)ij\ = 1/A^, so that C'i(p^) = N{N - 1 )T = 



° 2 1^1 4 0 2 1^1 4 


FIG. 4: Histogram of absolute values of rescaled 
off-diagonal elements y = pij of random mixed state. 
Left real state and right complex state. The solid (red) 
curves represent the distributions Xidl/I) (left) and 
V 2 x 2 {V 2 \y\) (right). 


N — 1. This result is greater, than the mean value (62) 


for a complex random state in a generic basis, by a factor 
7r/4 - 0.785. 

The random mixed state p = can be con¬ 

sidered as a normalized Gram matrix for an ensemble of 
N random complex vectors. In the asymptotic limit the 
normalization pre-factor A/TrGG^ tends to a constant. 
Thus in the limiting case we may consider 


1 ^ - 
Pi] = ^ GikGjk- 
^ k=l 


(63) 


From the central limit theorem we get that pij has a nor¬ 
mal distribution for N ^ 1. This implies that the ab¬ 
solute value \pij\ is described, up to a scaling factor, by 
an appropriate Xv distribution. It is the distribution of 
the square root of the sum of squares of independent ran¬ 
dom variables having a standard normal distribution m- 
The Xi/ distribution has a single parameter, the number 
v of degrees of freedom, equal to the number of summed 
normal variables. Hence, we get: 


Pmiy) = Xliy) = Y^e y'C 

Pciy) = V2x2{V2y) = 2xe-y\ 


(64) 


where y = y/Nlpijl- Fig. W shows that statistical distri¬ 
bution of modulus of off-diagonal elements of a random 
mixed state which can be asymptotically described by 
the X distribution. This property allows us to describe 
asymptotic behavior of the coherence of real and complex 
random mixed states. 


G?^(p“) - Vn xPudx = (65) 

Gf^ip^) Vn I xPcdx = VN^. ( 66 ) 


In the above, we have shown that the Li coherence, 
Gli , scales as the dimension N for generic pure states and 
as VN for typical mixed states. Note that expression ( 66 ) 
asymptotically holds also for contradiagonal states, for 
which Pa = 1/N. 
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IX. DYNAMICAL SYSTEM-COUPLED 
QUANTUM KICKED TOPS 

To show a direct link with the theory of quantized 
chaotic systems we analyze the model of quantum kicked 
top. It is described by the angular momentum operators 
Jx-,Jy,Jz-, satisfying standard commutation rules. The 
dynamics of a single system consists of a periodic unitary 
evolution in the Hilbert space of dimension N = 2j + 1, 
followed by an infinitesimal perturbation characterized 
by the kicking strength k. To analyze the effects of 
quantum entanglement the model consisting of two cou¬ 
pled tops described by the Hamiltonian, H{t) = Hi{t) + 
H- 2 (t) + Hi 2 {t), with Hi{t) = piJy, + 

Hi 2 {t) = e/j (g) JzJ and j = This 

model was studied in glHS]. 

We set the standard values of the parameters pi = 
P 2 = TT12 and consider the unitary one-step time evolu¬ 
tion operator, U = Ui 2 {Ui (g C/ 2 ), where 



Consider a pure separable state |0 C) |0 evolved unitarily 



1 


— r.' 

Dti 


D “ 

2 

D ^ 

^ _ 

0 

^ _ 


0 ^- 0 ^- 

0 50 i too 0 50 ^ 100 


FIG. 5: Trace distance Utt between initially orthogonal 
reduced states of the kicked tops as a function of time t 
obtained for ji = j 2 = 60. Convergence to the 
asymptotic value D represented by the horizontal line 
depends on a) kicking strength, k = 3.2 (upper dashed 
line), k = 3.5,3.8 and k = 4.0 (solid line) and b) the 
coupling constant, e = 0.005 (upper dashed line), 
e = 0.01, 0.1 and e = 1 (solid line). 


for t steps, which gives |'0() = U*{\1) 0 |Z)) for Z = 1,2. 
After the unitary dynamics we perform partial trace over 
the second subsystem, thus defining two states: 

ai{t) =Ti 2 \'tp\){'ip\\ and (72{t) = Ti 2 \i^\) (68) 


and study the spectrum of their difference F = cji — cr 2 - 
Level density of F obtained for fc = 6, corresponding 
to the chaotic regime, with Wigner level spacing statis¬ 
tics, e = 0.01 and four values of the dimensionality ratio, 
c = (2j2 -I- l)/(2ji -I- 1), can be described by distribu¬ 
tion (10) as shown in Fig. This observation confirms 
the conjecture that a quantized deterministic chaotic sys¬ 
tem may lead to generic, random states. Consider the 
trace distance UTr(cri(t), cr 2 (t)), which at t = 0 is equal 


to unity as the initial states are orthogonal. For larger 
times t the trace distance tends to the number D derived 
in ( [T6| ) and the convergence rate is faster for strongly 
chaotic systems (large k) and large coupling strength e - 
see Fig. 

The trace distance UTr(a’i(t),cr 2 (t)) between two ini¬ 
tially orthogonal states undergoing the dynamics of the 
coupled kicked top ( [^ decreases monotonously and 
tends to the universal number D. The rate of conver¬ 
gence to the asymptotic value depends on the parameters 
of the model, as shown in Fig. However, if we repeat 
the procedure taking for initial states the tensor products 
of two coherent states localized nearby, the trace distance 
between reduced states initially grows in time. This fea¬ 
ture is related to the fact that the dynamics taking place 
in the reduced system is not Markovian [44] and that 
there exist correlations between both systems. 

Note that the investigated state of the first subsystem 
is obtained by a reduction over the second subsystem, 
a{t) = TT2\'tp{t)}{'ip{t)\. Thus to obtain the next itera¬ 
tion a(t -I- 1), one needs to purify a{t) in a specific way 
to get the bipartite pure state \ip{t)), evolve it unitar¬ 
ily into -I- 1) = U\tp) and perform the partial trace, 
a{t + 1) = Tr2|V'(t + 1))(^(/ + 1)1- The purification pro¬ 
cedure is not unique, and the information encoded in a 
is not sufficient to recover the desired pure state |'0(t)). 
Hence the dynamics of the reduced state, a{t) a{t+l), 
is not Markovian. Furthermore, the quantum dynamics 
is explicitly defined for a full system consisting of two 
coupled tops, but the reduced dynamics of a single sub¬ 
system is not well defined independently of the correla¬ 
tions. 


X. CONCLUDING REMARKS 

We evaluated common distances between two generic 
quantum states. Due to the concentration of measure in 
high dimensions these distances converge to deterministic 
values - see Tab.jH] Our results are directly applicable for 
quantum hypothesis testing m as they imply concrete 
bounds on the distinguishability between generic states. 
Asymptotic value of the root fidelity, '/F = 3/4, can be 
used as a universal benchmark for future theoretical and 
experimental studies based on this quantity and Bures 
distance. 

Although obtained expectation values are exact in the 
asymptotic limit our numerical results presented in Fig.[^ 
show that they can be used for N of the order of ten. 
Furthermore, results obtained improve our intuition con¬ 
cerning the structure of the body Qjv of quantum states 
and describe coherence with respect to a generic basis 
of a typical pure and mixed quantum state. Moreover, 
we demonstrated that a typical mixed state of a large 
bipartite system is weakly entangled. 

Results presented in this paper are obtained with use 
of the free probability calculus and hold in the asymp¬ 
totic limit N —>■ 00 . It would be therefore interesting to 
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apply other methods to obtain finite size corrections to 
the expressions derived in this work and to show analyt¬ 
ically, for what system sizes the correction terms can be 
neglected in practice. In the case of the MV distribu¬ 
tion this has been studied in m, although these results 
do not directly apply in the case studied in this work. 
The free convolution works only in the asymptotic limit, 
thus obtaining a distribution similar to SMV for finite 
N requires careful treatment. The problem of estimat¬ 
ing the average coherence for quantum systems of a finite 
size was recently studied in [35] and their results are 
consistent with our asymptotic expressions. 


X 

Vx{p, ,vl) 

Dx(p,cr) 


Tr 

« 0.414 

47r 

K 4 « 0.568 

1 

HS 

0 

0 

1 

00 

0 

0 

1 

T 

Ti « 0.368 

1 

yih2 « 0.833 

B 

^20.550 

^ « 0.707 

V2 « 1.414 

E 

El « 0.518 

0.614 

yin2 « 0.833 

H 

^20.550 

^2-2(^)^«0.748 

« 1.414 


TABLE II: Typical distances between generic mixed 
states p and a and pure states \'ij)) and | 0 ) of a large 
dimension N. Definitions of D^r, Dt, Db and Dh are 
given in the text, while IIhs(p, ~ 

and Dao(p,a) = max^ \\i{p — cr)|, and the entropic 

distance [45| is D\{p,a) = H 2 ^ 5(1 — \/F{p, a))'j , with 

H 2 (x) = —xhix — (1 — a:) ln(l — x). Analytical formulas 
for the numbers Ti and Ei are provided in Appendix |C| 
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Appendix A: Limiting behavior of functions of two 
independent states 

In this section we consider the limiting value of the 
trace of a product of two functions on random states. 


Let p and cr be random mixed states generated accord¬ 
ing to the flat measure t'ns in Djv with eigenvalues de¬ 
scribed by the measures p!'^ = ^ ^N\i{p) and = 

For large N with probability one the mea¬ 
sures pn and vjq converge weakly to pmv ■ Because p and 
a are asymptotically free random matrices, thus for any 
analytic functions g and h the matrices g{p) and h{a) 
are also asymptotically free. The product g{p)h{a) has 
a limiting, non-random distribution of eigenvalues, given 
by an appropriate multiplicative free convolution. As a 
result the normalized trace of the product tends, almost 
surely, to a deterministic quantity, equal to the mean 

Trg(p)/i(cr) > m, (AI) 

Af—>00 

where m = limjv->.oo ^Trg{p)h{a). 

Now we will calculate the limiting value of the mean 
trace of the product. The eigendecomposition gives us 

ETrg{p)h{a)=EY, 9 Mp)]h[Xj {a)] (A2) 
ij 


where 17 is a unitary transition matrix from eigenbasis of 
p to eigenbasis of cr. The independence of eigenvectors 
distribution from eigenvalues distribution gives us 

ETrg(p)h(a) = E^ 5 [A(p)]h[A,(a)] (E|t7,,f) . (A3) 


Since the distribution of p and a is invariant with respect 
to unitary rotations, U has a Haar distribution, for which 
we have, by the permutation symmetry, E|t7y p = 
This gives us 

EFvg{p)h{a) = ^ETr 5 (p)ETrh(cT). (A4) 

Now we can make use of the weak convergence of the 
eigenvalues distributions similarly to ([^ and get 

ETrg(p)h(cr) > [ Bg^h{s,t)dpMv{t)pMv{s), 

N—^oo I 

(A5) 

where 

i?,,,(s,f)=Jim^iVg(^)h(^). (A 6 ) 


Appendix B: Distances as a function of the 
rectangularity parameter 

In this section we present the average trace and the 
rescaled Hilbert-Schmidt distance between random states 
distributed with respect to the induced measure vk- It 
is convenient to parameterize the measure by the rect¬ 
angularity parameter c = K/N, where N stands for the 
dimension of the principal system, while K denotes the 
dimension of the environment. 
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(a) Dxr 



(b) VNDhs 


FIG. 6: Trace distance (a) and rescaled Hilbert-Schmidt 
distance M as a function of the rectangularity 
parameter c. So id blue curve denotes the distance of a 
random state p from the maximally mixed state p*. 
Dashed red curve shows the distance between two 
random states. The vertical lines marks c = 1. The 
solid horizontal line marks D, the distance between two 
generic states for c = 1. The dashed horizontal line 
marks a, the distance of a generic state form the 
maximally mixed state p*. 


Following (§ we integrate the function \t — c\ over 
the Marchenko-Pastur distribution AiVc to get the av¬ 
erage trace distance IlTr(p, p*) between a random state 
p and the center of the set, p* = 1/N. In analogy 
to Eq. (161 computing the average |a;| with respect to 
the symmetrized MP distribution SMVc one obtains the 
mean value of the trace distance IlTr(p, cr) between two 
random states. 

These quantities for the trace distance are shown in 


in Fig. 6a as a function of the rectangularity parameter 
c. In the limit c —> oo the distribution SMVc tends 
to the circular law rescaled by so the distance 

DTr(p, cr) tends to zero as 4,^/2 j{iiry/c). Note that above 
average improves the bound recently derived [46] in order 
to show that exponential decay of correlations implies 
area law. Integrating the function \t — c| over the MVc 
distribution we obtain that the average trace distance 
from the maximally mixed state p*, D'Yrip^p*) tends to 
zero as 4/{Z■K^/c). 

In order to calculate the rescaled Hilbert-Schmidt dis¬ 
tance between two generic states, y/N Dhs (Pj cr); we eval¬ 
uate the second moment of the SMVc distribution. In 
the limit c —>■ oo the distance tends to zero as \f2ly/c. 
To obtain average Hilbert-Schmidt distance of a generic 
state from the maximally mixed state, •\/]VI?hs(P) P*\ we 
integrate the function — over the Marchenko-Pastur 
distribution MVc- In t he l imit c —>■ oo this distance be¬ 
haves as lly/c- see Fig. 6b Note that for both distances 
in the limit c — >■ oo we get Dx{p, a) = \/2D^{p, p*), where 
X stands either for the trace or the Hilbert-Schmidt dis¬ 
tance. 


Appendix C: Transmission distance and entropic 
distance 


We provide here an analytical expression for two dis¬ 
tances listed in Tab. Square root of quantum Jensen- 
Shannon divergence (QJSD) leads to the transmission 
distance [12], I^t(p, o’) = ^QJSD(p, cr). Taking one ran¬ 
dom state p and the maximally mixed state p* we find 
that for large dimension N this distance converges to 


Ti-[DTip,P*)f^l + ^ + ^log2+ 


log '^4870847 - 2178309v^, 


so that Ti ~ 0.368. In a similar way we find the entropic 
distance defined in the caption to the Table I, 


El = De{p,P*) 


(C2) 


\- 


STrlog - 16coth ^ 




6tt 


= 0.518. 
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